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Abstract 

We have observed the transversal vibration mode of suspended carbon nanotubes at millikelvin tem- 
peratures by measuring the single-electron tunneling current. The suspended nanotubes are actuated 
contact-free by the radio frequency electric field of a nearby antenna; the mechanical resonance is de- 
tected in the time-averaged current through the nanotube. Sharp, gate-tuneable resonances due to the 
bending mode of the nanotube are observed, combining resonance frequencies of up to pq = 350 MHz 
with quality factors above Q = 10 5 , much higher than previously reported results on suspended carbon 
nanotube resonators. The measured magnitude and temperature dependence of the Q-factor shows a 
remarkable agreement with the intrinsic damping predicted for a suspended carbon nanotube. By ad- 
justing the RF power on the antenna, we find that the nanotube resonator can easily be driven into the 
non-linear regime. 

High-quality resonating systems, providing high frequency resolution and long energy storage time, play 
an important role in many fields of physics. In particular in the field of nanoelectromechanical systems [1, 2], 
recent research has led to the development of high-frequency top-down fabricated mechanical resonators 
with high quality factors [3, 4, 5, 6]. However, when miniaturizing mechanical resonators to make them 
lighter and to increase their resonance frequency [1], the quality factor tends to decreases significantly from 
surface effects [2]. High Q- values combined with high resonance frequencies are an important prerequisite 
for applications such as single- atom mass sensing [7, 8, 9] and fundamental studies of the quantum limit 
of mechanical motion [10]. Single- wall carbon nanotubes present a potentially defect-free nanomechanical 
system with extraordinary mechanical properties: in particular the high Young's modulus (E = 1.2 TPa) in 
combination with a very low mass density (p — 1350 kg/m 3 ) [8, 11, 12, 13]. While these favorable properties 
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should result in quality factors of the order of 2 x 10 5 [14], the observed Q-factors of nanotube resonators both 
at room temperature [11, 12, 15, 16] and in low temperature experiments [7, 8] have not exceeded Q <~ 2000. 

Here we report on the observation of mechanical resonances of a driven suspended carbon nanotube at low 
temperatures with quality factors above 10 5 and resonance frequencies ranging from 120 MHz to 360 MHz. 
The resonances are detected with a novel detection scheme which uses the non-linear gate-dependence of the 
current through the suspended nanotube quantum dot. In addition, we show that the nanotube resonator 
can easily be tuned to the non-linear regime, and that the operating temperature affects the non-linearity 
and the quality factor of the resonator. 

Suspended carbon nanotube devices are made by growing nanotubes between platinum electrodes over 
an 800 nm wide pre-defined trench. The device geometry is shown in Figure 1(a). The fabrication method is 
discussed in detail by Steele et al. [17]. There, the device includes three local gates for tuning the confinement: 
here, however, we apply the same voltage to all three gates, so that they act together as one single gate. Since 
no device processing takes place after nanotube growth and the entire device is suspended, the nanotubes 
arc highly defect-free and do not suffer from potential irregularities induced by the substrate surface [17, 18]. 
The fabrication method also offers the advantage that the resonator is not contaminated with resist residues. 

After fabrication, the suspended nanotube devices are mounted in a dilution refrigerator with filtered 
twisted pair cabling attached to source, drain, and gate contacts (see Figure 1(a)). This configuration allows 
us to apply dc gate and bias voltages to the suspended nanotube, and measure the current flowing through it. 
To minimize heating, we drive the nanotube resonator with the electric field radiated from a radio frequency 
(RF) antenna positioned near the sample (<~ 2 cm) instead of connecting high-frequency cables directly 
to the sample. Measurements arc performed at temperatures down to the base temperature of the mixing 
chamber of the dilution refrigerator, Xmc — 20 mK. 

Figure 1(b) shows the Coulomb oscillations of a semiconducting carbon nanotube with a suspended length 
of 800 nm. A highly regular addition spectrum with clear four-fold degeneracy is visible, characteristic for 
a defect-free single- wall carbon nanotube [19, 20]. From the magnetic-field dependence of the position of 
the Coulomb oscillations close to the semiconducting gap [21], we find the radius of the nanotube, r, to be 
between 1-1.5 nm. The value of the semiconducting gap « 0.3 eV is estimated from the gate range between 
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electron and hole conduction in the device at low temperatures, which is in agreement with this value for r. 

When an ac voltage Vr F with frequency v is applied to the antenna, we observe a resonant feature at 
a well-defined frequency in the dc current flowing through the nanotube. Figure 1(c) shows an example 
of such a measurement at a large radio frequency (RF) voltage, or equivalently a high generator power. A 
sharp resonant feature is clearly visible at v = 294 MHz. Zooming in on this feature at a lower power (Figure 
1(d)) reveals a resonance peak with a narrow lineshape. A numerical fit of this data yields a quality factor 
Q = 140670 (see below for a discussion of the expected lineshape) . We have also performed measurements 
on a second device displaying similar resonant peaks with Q-factors up to 20000; the results on that device 
are shown in the Suppl. Information. 

The resonance observed in parts (c) and (d) of Figure 1 can be attributed to the flexural vibration mode 
of the suspended nanotube [8, 7, 11, 12]. To verify this, we electrostatically induce tension in the nanotube 
by applying a dc gate voltage V g to the back gate electrode [11, 12]. The dc gate voltage dependence is 
shown in Figure 2. When decreasing the gate voltage from zero to more and more negative values, the 
resonance is tuned to higher frequencies by almost a factor of three: from less than v a = 140 MHz at 
V g = — 1 V to vo = 355.5 MHz at V g = —6.5 V. For the latter resonance frequency, the thermal occupation 
[10] n = 1/2+ [expihvo/knTMc) — I] -1 would be 1.2 at 20 mK, suggesting that the resonator would be close 
to its quantum ground state in the absence of the driving fields required for our detection scheme. 

We have extracted the resonance peak positions from the data in Figure 2 and plotted them in the 
inset. The red line shows the gate dependence of the resonance frequency calculated with a continuum 
model for the fundamental flexural bending mode [22, 12, 23]. The parameters are ^bonding = 132.0 MHz, 
V* = -2.26V, and T = 0, where ^bending is the resonance frequency in absence of residual tension To and 
Vg* marks the cross over between the weak and strong bending regime. At high gate- voltage the model 
calculation deviates slightly from the experimental values. This is so far not fully understood and may be 
related to large static displacements of the nanotube in a complex electrostatic environment [24]. The value 
^bonding = 22.4/27rL 2 • r^E/p = 132.0 MHz, assuming a tube length L = 800 nm, yields a nanotube radius 
of 1.6 nm, in good agreement with the band-gap and magnetic field estimates. 

Depending on the gate voltage, the resonance either appears as a dip (Figure 3(a)) or as a peak (Figure 
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3(b)). Dips are found around the maxima of the Coulomb oscillations; away from these maxima peaks are 
observed. This indicates that the detection of the mechanical modes is due to electrostatic interactions as 
we will now show. We model the effect of a small change in gate voltage SV g on the current flowing through 
the nanotube by a Taylor expansion of I(V g + SV g ) around SV g — 0. A crucial point in this expansion is that 
the second (and third) order term cannot be neglected, since the current flowing through the nanotube is 
strongly non- linear in the vicinity of the Coulomb oscillations. This is in contrast to the mixing technique 
[11, 12], where only the linear term in the expansion is needed. 

The motion of the nanotube enters the measured current as follows: On resonance, the nanotube position 
u(t) = u cos {2-Kv t) oscillates with a finite amplitude uo, which periodically modulates the gate capacitance 
C g by an amount C g c = (dC g /du) u . The current flowing through the nanotube does not just depend on the 
gate voltage itself; more specifically, it depends on the product of the gate voltage and the gate capacitance 
- the so-called gate-induced charge [11, 12, 25]. A modulation of the capacitance due to the motion of the 
nanotube therefore has the same effect on the current as if an effective ac gate voltage Vg ac,cff = V g C g c /C g 
were applied to the gate-electrode. The time-dependent current can then be calculated by inserting SV g = 
y g ac ' eff cos(27r^t) into the Taylor expansion of I(V g + 5V g ). 

Since the mechanical resonance frequency is much larger than the measurement bandwidth, time-averaged 
currents are detected in our setup. We find that the time-averaged mechanically induced current equals: 



where only even powers of u enter the low-frequency current due to averaging. The change in dc current 
on mechanical resonance AI = I — I is thus proportional to the local curvature d 2 I/dV g of the Coulomb 
blockade oscillations I(V g ). 

Using measured Coulomb oscillation traces where no driving signal was applied (black line in Figure 3(c)), 
we have numerically calculated the behavior of a current time-averaged over V g c ' cS . The result is shown as a 
red line in Figure 3(c). Figure 3(d) shows the difference AI between the time-averaged and the static current 
of Figure 3(c). On top of a Coulomb oscillation, the curvature is negative and the averaged current (red) 
is smaller than the static current (black), resulting in a dip in the current on resonance, when mechanical 
motion takes place. On the other hand, AI is positive on the flanks of the Coulomb oscillations as the 
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curvature is positive there. This can be compared with the traces I{v) shown in Figure 3(a) and (b), and 
with the measurements of AI shown in Figure 3(e). Here we plot the amplitude AI(i/ ) of the mechanical 
response in the dc current I{u) for different gate voltages. The gate voltage dependence of the extracted 
amplitude values in dc current is in good qualitative agreement with the predictions of the model as shown 
in Figure 3(d). 

The model also allows for a quantitative analysis of the peak shape and for an estimate of the displacement 
amplitude u in the case of resonant driving, by evaluating the change in dc current AI alone. We first note 
that u can be described by the response of a damped driven harmonic oscillator [12, 22]. From Eq. 1, we 
see that AI cx u 2 , so that the measured mechanical response (dip or peak) in the current is given by the 
square of the harmonic oscillator response function (SHO). For the resonance presented in Figure 1(d), we 
find v — 293.428 MHz and Q = 140670. This Q-value is nearly two orders of magnitude higher than previous 
reported values of the flcxural vibration modes in nanotubes [7, 8, 11, 12]. Such high Q- values make this 
type of device very suitable for mass detection. From the measured response in Figure 4(e) we estimate (see 
Suppl. Information) a mass sensitivity of 7yg/v / Hz, i.e., in one second it should be possible to determine if, 
for example, a He atom has adsorbed onto the nanotube. 

The displacement amplitude u in the case of resonant driving is estimated by modelling the capacitance 
between the nanotube and the back-gate as the capacitance between an infinite wire and an infinite conductive 
plane [22]. Using a device length of L — 800 nm, a tube radius r =1.5 nm and a gate distance h = 230 
nm, we obtain C g = 7.8 aF and dC g /du = — 5.9zF/nm. The calculated capacitance value is consistent with 
the experimentally determined value of C g = 8.9 aF as determined from the Coulomb peak spacing. For the 
resonance in Figure 4(b), with d 2 I/dV 2 = 4A3^A/V 2 and AI(i/ ) =1.05 pA, we estimate the oscillation 
amplitude of the nanotube to be «o(^o) = 0.25 nm on resonance. This amplitude is two orders of magnitude 
larger than that of the thermal fluctuations {k^,T/2 = m{2-KU )' 2 u 2 h /2) of the nanotube [1], which is ~ 6.5 pm 
at 80 mK, and its estimated zero-point motion [10, 26] of 1.9 pm at this gate voltage. 

When driving the nanotube resonator with large antenna voltages, we consistently observe hysteretic peak 
shapes and a strong frequency pulling of the resonance peaks (i.e. the frequency decreases for a larger motion 
amplitude [27, 28]). Figure 4(a)-(d) show examples of the shape of the resonance peak at V g — —5.16V and 
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V s d = 0.35 mV for four different driving powers. Black lines indicate the sweep direction with increasing 
frequency; red lines the one with decreasing frequency. At the lowest power, the mechanical resonance peak 
is not visible in the noise. With increasing driving power the resonance peak first shows a linear response 
with its characteristic SHO shape (Figure 4(b)). At higher powers hysteresis sets in, which becomes more 
pronounced with increasing RF power. This bistability is consistent with what is expected for a non-linear 
mechanical (Duffing) resonator [1, 27]. 

We have studied the dynamic range [29, 2, 24, 30] in more detail and found that the driving powers where 
the (linear) peak disappears in the noise and where nonlincarity sets in depend on the temperature. An 
example of this effect is shown in Figure 4(e)-(h). These panels show that for a fixed gate voltage and driving 
power, the nanotube resonator response changes from non- linear to linear when the operating temperature 
is increased from 20 mK to 160 mK. This temperature-dependent behavior hints at a decrease in Q-factor as 
the temperature is increased. 

To study the temperature dependence of the quality factor in more detail, we have determined Q at 
different temperatures. For a gate voltage of -5.16 V, three examples of resonance traces are depicted in 
Figure 5(a)-(c). Note that because the dynamic range is temperature dependent, the RF power is adjusted 
at every temperature to ensure a linear response. In Figure 5(d), we plot the Q-factor extracted in the 
linear regime for eight different temperatures in the range 20 mK< T MC < 1 K. The error margins are 
estimated from ensembles of responses at the same temperature. The Q-factor changes by a factor four in this 
temperature range. At the lowest temperatures, the Q-factor reproducibly reaches values above 10 5 . These 
lowest temperature values are close to the intrinsic Q-values calculated with molecular dynamics simulations 
on single- walled carbon nanotube oscillators [14]. Interestingly, these calculations predict a T~°- 36 power law 
dependence of the Q-factor with temperature. The red line in Figure 5(d) shows this dependence; the data 
is consistent with this prediction. This T~ 36 dependence has also been observed at low temperatures in 
top-down fabricated devices [31, 6]. Note that the Q-values of our nanotube resonator are much higher than 
the ones following the trend of the volume surface ratio in top-down fabricated devices [2]. 

In conclusion, using a novel detection mechanism, we have measured the bending mode resonance of 
suspended carbon nanotubes in the single-electron tunneling regime. Sharp gate-tunable resonances are 
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found with high Q- values (Q > 10 5 ), which can easily be driven into the nonlinear regime by increasing the 
driving power on the RF antenna. By inducing tension with a gate voltage the frequency can be tuned above 
350 MHz, so that the thermal occupation of the resonator approaches 1 . Shorter devices should have even 
higher resonance frequencies corresponding to temperatures higher than the mixing chamber temperature of 
the dilution refrigerator. These resonators are therefore in their quantum mechanical ground state, which 
opens up the way to new exciting experiments on the quantum aspects of mechanical motion. 
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Figure 1: (a) Schematic drawing of the chip geometry, antenna, and measurement electronics. The nanotube 
acts as a doubly clamped beam resonator, driven by an electric field E(t). The displacement of the nanotube 
is u(t). (b) Example trace of the dc current at V^d = 50 /iV as a function of gate voltage, demonstrating the 
regularity of the Coulomb peaks. It shows the four-fold degeneracy typical for clean single-wall carbon nanotubes. 

(c) When the frequency v of an RF signal on the antenna is swept with fixed V g and V s d, a resonant peak emerges 
in I(v). An example of such a resonance is shown for a driving power of — 17.8dBm at a temperature of 20 mK. 

(d) Zoom of the resonance of (c) at low power (— 64.5dBm). The red line is a fit of a squared damped driven 
harmonic oscillator response to the resonance peak. For both (c) and (d) V g = —5.16 V and V s d — 0.35 mV. 
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Figure 2: \dl /dv\ as a function of frequency v of the ac voltage on the antenna and the dc gate voltage V g on 
the back-gate electrode. Horizontal stripes are caused by electrical (cable) resonances [11, 12]; the narrow vertical 
stripe pattern is related to the Coulomb blockade oscillations. In addition, a gate-dependent resonant feature is 
clearly visible. Inset: Comparison of the extracted resonance frequency to the continuum model for the bending 
mode with lending = 132.0 MHz, V* = 2.26V, T = 0, and a shift of 0.775V in gate voltage to account for 
an offset in the charge neutrality point of the nanotube from V g = 0V and the band gap region [12, 22]. The 
parameters are discussed in the text. An apparent shift of the mechanical resonance frequency at Uq ~ 230 MHz 
is caused by an electrical (cable) transmission resonance, leading to a strong increase in transmitted RF power 
and distorted peak shapes. 
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Figure 3: Averaging model for the current at resonance, (a), (b) Measured frequency sweeps demonstrating 
the sign change of the resonance amplitude depending on the gate voltage (RF power — 13dBm, V^d = 0.1 mV, 



-5.17V (a) and V g = —5.16V (b)). (c) The black line shows the measured dc current as function of 



gate voltage I(V g ) for V a d = 0.1 mV (no RF). The red line, shows the effect of an (effective) ac gate voltage on 
the dc current. This average current (Eq. (1)) is calculated using the measured data and V g ac ' eff = 2mV. (d) 
Predicted resonance signal amplitude AI calculated by subtracting the dc current from the current averaged over 
an effective gate voltage Vg ac ' cff = lmV. For a small Vg C ' eS , the signal is proportional to the second derivative 
d 2 I/dV„ of the black trace shown in (c), as expressed in Eq. (1). At the top of the Coulomb peak, AI is negative, 
whereas on the flanks of the Coulomb peak, it is positive. Note that in (c) a larger value of V g c ' was used to 
exaggerate the difference between the black and red curves for illustrative purposes, (e) The measured resonance 
peak amplitudes obtained from I(v) traces similar to Figure 1(d), for V sc i = 0.1 mV and a RF power of —48 dBm. 
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Figure 4: Evolution of the resonance peak with increasing driving power (a-d) and temperature (e-h). Black (red) 
traces are upward (downward) frequency sweeps, (a) At low powers, the peak is not visible, (b) Upon increasing 
power, a resonance peak with Q=128627 appears. (c,d) As the power is increased further, the lineshape of 
the resonance takes on a non-linear oscillator form, with a long high frequency tail and a sharp edge at lower 
frequencies. It also exhibits hysteresis between the upward and downward sweep that increases with driving power, 
characteristic of a non-linear oscillator. The traces (a-d) are taken at 80 mK. (e-h) Forward (black) and reverse 
(red) frequency sweeps at a fixed driving power as a function of temperature. At low temperatures, the peak 
shape is non-linear and strongly hysteretic. At the same power, but higher temperature, the amount of hysteresis 
decreases significantly. At a temperature of 160 mK, hysteresis and asymmetry are no longer apparent; at the 
same time, the signal amplitude (and with it also signal to noise ratio) is decreased, suggesting a decrease in the 
Q-factor with increasing temperature. The working point of traces (a-h) is at V g — —5.16 V and V^d = 0.35 mV. 
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Figure 5: Temperature dependence of the Q-factor. (a-c) Fits of a squared harmonic oscillator response to the 
resonance in the linear regime at low powers for different temperatures at V g — —5.16 V and V s d — 0.35 mV. 
(d) A plot of the Q-factor vs. temperature obtained from linear response traces. Q decreases with increasing 
temperature. The red line shows a T~ - 36 power law dependence (see text). 
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1 Mass sensitivity 

The mass sensitivity is estimated from the data plotted in Figure 4(e) of the main text. In the experiments, 
the mass sensitivity is limited by the current noise, which has a spectral density S^ 2 = 0.12pA/-\/Hz. The 

1 /2 

mass sensitivity Sm can then be calculated as follows: an added mass 5m on the nanotube changes the 

resonance frequency by: 

5vo = w^-S™ = 7r-Sm, (SI) 
g dm 2m 

i-^j where m = 5.1 x 1CP 21 kg is the mass of a 800 nm long single-walled nanotube with a 1.5 nm radius. When 

S3 

the resonance frequency shifts, the current through the nanotube is modified by: 

(N 5I=^5v ~-^-8v Q . (S2) 

CO 

The latter approximation, which is valid for a high Q resonator, allows us to relate the change in current to 
the measured slope of the response function. For the data in Fig. 4(c) the slope of the red line, just right of 
the jump is dl/dv = 6.0 x 10~ 16 A/Hz. The mass sensitivity is then calculated from sU 2 = §^ (f^)" 1 S 1 / 2 , 

o 1/2 

L; which yields = 7.0yg/VHz = 4.2u/VHz. Here, u is the (unified) atomic mass unit, so it is possible to 

detect a mass change as small as a single helium atom within one second. 



*Present address: Institute for Experimental and Applied Physics, University of Rcgensburg, 93040 Regensburg, Germany 
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2 Device B 



Figure 1 and Figure 2 show measurements on a second device. This device also has a suspended length 
of 800 nm. From room temperature measurements it is inferred that device B is a large (E g /kb > 300 K) 
bandgap nanotube as well. 




263.8 263.9 263.8 263.85 
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Figure 1: Examples of measured resonances in device B in the linear (a) and non-linear regime (b) at 20 mK. 
Settings: V g = -4.241V, V sd = 2.0 mV, RF power -47dBm in (a) and V g = -4.241V, V sd = 1.5 mV, RF 
power — 44.5dBm in (b). 
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Figure 2: \dl /dv\ as a function of frequency v of the ac voltage on the antenna and the dc gate voltage V s 
on the back-gate electrode for device B. Left of the dashed line a source-drain voltage of 4mV was used; on the 
right side V s d = lOmV. The RF power was — 13dBm everywhere. Inset: Comparison of the extracted resonance 
frequency to the continuum model for the bending mode with ^bending = 193.7MHz, V* — 4.14 V, T = and 
a horizontal offset of 1.65 V to account for a shift in the charge neutrality point and the band gap region. 
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